Pairing state in the rutheno-cuprate superconductor RuSr 2 GdCu20 8 : 
A point contact Andreev reflection spectroscopy study 



S. Piano, F. Bobba, F. Giubileo and A. M. Cucolo, 
M. Gombos, A. Vecchione 

Physics Department and INFM-CNR SUPERMAT Laboratory, 
University of Salerno, Via S. Allende, 84081 Baronissi (SA), Italy 
(Dated: February 2, 2008) 

The results of Point Contact Andreev Reflection Spectroscopy on polycrystalline RuSr2GdCu2 08 pellets 
are presented. The wide variety of the measured spectra are all explained in terms of a modified BTK model 
considering a d-wave symmetry of the superconducting order parameter. Remarkably low values of the energy 
gap A = (2.8 ± 0.2)meV and of the 2A/fc_BT c ~ 2 ratio are inferred. From the temperature evolution of the 
dl /dV vs V characteristics we extract a sublinear temperature dependence of the superconducting energy gap. 
The magnetic field dependence of the conductance spectra at low temperatures is also reported. From the A vs 
H evolution, a critical magnetic field H C2 ~ 30T is inferred. To properly explain the curves showing gap-like 
features at higher voltages, we consider the formation of a Josephson junction in series with the Point Contact 
junction, as a consequence of the granularity of the sample. 
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I. INTRODUCTION 

Point Contact Spectroscopy is a versatile technique widely 
used to study the basic properties of superconductors, such 
as the density of states at the Fermi level and the supercon- 
ducting energy gap. The technique consists in establishing a 
contact between a tip of a normal metal (N) and a supercon- 
ducting sample (S), thus forming a small contact area that is a 
"Point Contact" junction. By varying the distance and/or the 
pressure between tip and sample it is possible to obtain dif- 
ferent tunnel barriers, that is different conductance regimes. 
Indeed, quasiparticle tunnel spectroscopy is obtained for high 
barriers, while Point Contact Andreev Reflection (PCAR) 
spectroscopy is achieved in case of low barriers. Often in 
the experiments, intermediate regimes are realized, in which 
through the N/S contact both quasi-particle tunneling and An- 
dreev reflection processes occur. 

Andreev Reflections 2 - 3 take place at the N/S interface when 
an electron, propagating in the normal metal with an energy 
lower than the superconducting energy gap, enters in the su- 
perconductor forming an electron pair (Cooper pair) while a 
hole, with opposite momentum with respect to the incident 
electron, is reflected in the normal metal. A single reflection 
corresponds to a net charge transfer of 2e, where e is the elec- 
tron charge, from the normal metal to the superconductor. In 
the limit of low barriers at low temperatures, all the incident 
electrons at the N/S interface with energy eV < A are An- 
dreev reflected and the conductance doubles the normal states 
value. 

In this paper we report on PCAR studies carried out in the 
hybrid rutheno-cuprate RuS^GdC^Os (Ru-1212) system 4 . 
This compound has recently drawn great attention among the- 
orists and experimentalists in the field of solid state physics 
due to the coexistence at low temperatures of superconduct- 
ing and magnetic ordering 5 . The Ru-1212 structure is similar 
to that of YBa2Qi307 with magnetic (2D) RuC>2 planes sub- 
stituting the (ID) Cu-O chains. The superconducting critical 
temperature in this compound strongly depends on the prepa- 



ration conditions with some reports showing transition onset 
as high as 50K 6 . The Ru-1212 also shows a magnetic phase 
below 135K. It has been reported that the magnetic order of 
the Ru moments is predominantly antiferromagnetic along the 
c axis 7 , while a ferromagnetic component has been observed 
in the RuC>2 planes, that act as charge reservoir 8 . At the mo- 
ment, due to complexity of this compound, an exhaustive de- 
scription of the interaction between the magnetic and super- 
conducting layers is still missing as well as an unambiguous 
evaluation of the symmetry of the energy gap. 

The paper is organized as follows: in Sec. [H] we briefly re- 
view the results of the BTK model 9 for a conventional s-wave 
superconductor and the recent extensions 10 for anisotropic s- 
wave and d-wave symmetry of the order parameter. In Sec.lIIII 
we describe the point contact experiments in polycrystalline 
Ru-1212 pellets showing a variety of conductance curves ob- 
tained at T = 4.2K. Satisfactory theoretical fittings are 
achieved by using a modified BTK model for d-wave sym- 
metry of the order parameter. In Sec. II VI we show that, due to 
the granularity of the samples, in some cases, the formation of 
a Josephson junction in series with the N/S contact occurs. In- 
deed, conductance curves showing gap-like features at higher 
voltages and dips in the spectra are well explained by this as- 
sumption. In Sec.|V]we report the temperature evolution of the 
conductance curves of a very stable junction. All the spectra 
are well reproduced by the d-wave modified BTK model, and 
we infer from the experiments the temperature dependence of 
the superconducting energy gap A. In Sec. lVII we analyze the 
magnetic field behavior of the measured conductance spec- 
tra, providing an estimation of the upper critical field of the 
Ru-1212 compound. Finally, in Sec. IVIII we summarize our 
results and draw some conclusions. 



II. THE BTK MODEL AND ITS EXTENSION 

In this section, for sake of clearness, we review the main 
results of the original Blonder- Tinkham-Klapwijk (BTK) the- 
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oretical model 9 , as developed for electronic transport through 
a Point-Contact junction between a normal metal and a con- 
ventional BCS superconductor. We also summarize the 
Kashiwaya-Tanaka 10 extension for asymmetric s-wave and d- 
wave superconductors. Indeed, a close comparison of the cal- 
culated conductance spectra is useful for a better understand- 
ing of the peculiar transport processes that occur at an N/S 
interface depending on the symmetry of the superconducting 
order parameter. 

Following the original paper, we write the expression of the 
differential conductance characteristics for a N/S contact that, 
according to the BTK model 9 , is given by: 
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where eV is the applied potential, Gnn = 4/(4 + Z 2 ) is 
the normal conductance expressed in term of Z, a dimen- 
sionless parameter modeling the barrier strength, f(E) is the 
Fermi function and A(E) and B(E) are, respectively, the An- 
dreev reflection and normal reflection probabilities for an 
electron approaching the N/S interface. Eq.Q shows that, 
while ordinary reflections reduce the transport current through 
the junction, Andreev reflections increase this by transferring 
two electrons (Cooper pair) in the superconducting electrode 
on the other side of the barrier. The case Z = corre- 
sponds to a completely transparent barrier so that the trans- 
port current is predominantly due to Andreev reflections and 
G NS (V < A)/G NN (V » A) = 2 is found (Fig[TJi). By 
increasing Z, the Andreev reflections are partially suppressed 
and the conductance spectra tend to the case of a N/I/S tunnel 
junction showing peaks at eV — ±A (Figs.^i,c). 

Recently, Kashiwaya and Tanaka 10 extended the BTK 
model by considering different symmetries of the order pa- 
rameter. Indeed, for a d-wave superconductor, the electron- 
like and hole-like quasiparticles, incident at the N/S interface, 
experience different signs of the order parameter, with for- 
mation of Andreev Bound States at the Fermi level along the 
nodal directions. The presence of Andreev Bound States mod- 
ify the transport current and the expression of the differential 
conductance is given by: 
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FIG. 1: Conductance characteristics, at low temperatures, for dif- 
ferent barriers Z as obtained by the BTK model for a Point Contact 
junction between a normal metal and a s-wave (a,b,c), a d-wave(d,e,f) 
and an anisotropic s-wave superconductor (g,h,i). 
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So, at a given energy E, the transport current depends both 
on the incident angle ip of the electrons at the N/S interface as 
well as on the orientation angle a, that is the angle between 
the a-axis of the superconducting order parameter and the x- 
axis. When applying Eqs. l|2j-(|6} to PCAR experiments, there 
is no preferential direction of the quasiparticle injection angle 
<p into the superconductor, so the transport current results by 
integration over all directions inside a semisphere weighted 
by the scattering probability term in the current expression. 
Moreover, because our experiments deal with policrystalline 
samples, the angle a is a pure average fitting parameter, which 
depends on the experimental configuration. 

In case of d-wave symmetry, for Z — > 0, the conductance 
curves at low temperatures show a triangular structure cen- 
tered at eV = 0, quite insensitive to variations of a with 
maximum amplitude Gns{V — Q)/Gnn(V >> A) = 2 
(Fig. ^1). However, for higher barriers, the conductance char- 
acteristics show dramatic changes as function of a. In par- 
ticular as soon as a ^ 0, the presence of Andreev Bound 
States at the Fermi level produces strong effects more evident 
along the nodal direction {a = 7r/4) for which Gns(V = 
0)/G NN (V » A) > 2 is found (Figs.[]>,f). 

For comparison, we report the conductance behavior for 
anisotropic s-wave superconductor, in which only the ampli- 
tude of the order parameter varies in the fc-space, while the 
phase remains constant and Eq. reduces to: 



Acos[2(a - ip)]. 



(7) 
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Again, in the limit Z — > 0, an increase of the conductance 
for E < A with a triangular profile is found with maximum 
amplitude G NS (V = 0)/G NN (V » A) = 2 at zero bias 
(FigH)- On the other hand, for higher Z, we obtain tunneling 
conductance spectra that show the characteristic "V"-shaped 
profile in comparison to the classical "U"-shaped structure 
found for an isotropic s-wave order parameter (Figs^i,i). We 
notice that in this case all the curves are quite insensitive to 
variation of the a parameter and a zero bias peak is obtained 
only for low barriers. 



III. PCAR SPECTROSCOPY ON RUSR 2 GDCU 2 8 : 
EXPERIMENTS AND THEORETICAL FITTINGS 

The Ru-1212 samples used for this study were direction- 
ally solidified pellets, grown by means of the Top-Seeded 
Melt-Textured method starting from Ru-1212 and Ru-1210 
(RuSi^GdOe) powder mixtures with a ratio Ru-1212/Ru-1210 
= 0.2. The details of the preparation procedure are reported 
elsewhere^ In the X-Ray diffraction patterns, a single Ru- 
1212 phase was found. In the resistivity measurements ver- 
sus temperature, the onset of the superconducting transition 
was observed at T°™ ~ 43X with T c (p = 0) ~ 2AK and 
AT C = 12K (AT C is defined as the difference between the 
temperatures measured at 90% and 10% of the normal state 
resistance). We notice that a broadening of the superconduct- 
ing transition is often observed in polycrystalline samples and 
it is usually related to the formation of intergrain weak Joseph- 



son junction 
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We address this point in the next section. 



To realize our experiments we used a Pt-Ir tip, chemically 
etched in a 40% solution of HC1, while Ru-1212 samples were 
cleaned in an ultrasound bath in ethyl alcohol. Sample and 
tip were introduced in the PCAR probe, in which three mi- 
crometric screws are allocated, each driven by its own crank. 
Two screws allow to vary the distance between tip and sam- 
ple, with a precision of 1/im and 0.1pm, respectively. The 
third screw is devoted to change the inclination of the sample 
holder varying the contact area on the sample surface. The 
Point Contact junctions were formed by pushing the Pt-Ir tip 
on the Ru-1212 pellet surface with the probe thermalized in 
the liquid He 4 bath. The Current- Voltage (/ vs V) charac- 
teristics were measured by using a conventional four-probe 
method and a lock-in technique with an amplitude of the ac 
current less than lpA was used to measure the differential 
conductance (dl/dV vs V) spectra as function of the applied 
voltage. 

In Fig. |3 we show a variety of normalized conductance 
spectra obtained at T — A.2K by establishing different con- 
tacts on different areas of the same Ru-1212 pellet. The junc- 
tion resistances varied between 10f2 and 100ft. By using the 
Sharvin relation, 14 it has been possible to achive an estima- 
tion of the size of the contact area. Indeed R = pi /4a 2 , 
where p = QAmVt cm is the low temperatures resistivity and 
I « 1000 A, as estimated in Ref. 1 1 . In our case, we have found 
that the typical contact size varied between 300A and 1000A. 

We observe that all the reported spectra are characterized 
by a Zero Bias Conductance Peak (ZBCP) with a triangular 
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FIG. 2: (Color online) The dl /dV vs V characteristics measured 
in different Ru-1212/Pt-Ir PC junctions at 4.2K. The experimental 
data (dots) are shown together with the best theoretical fittings (solid 
lines) obtained by a modified BTK model for a d-wave symmetry of 
the superconducting order parameter. 



structure, the main features appearing for each contact with 
different shapes, amplitudes and energy widths. Quite often, 
oscillations are observed on the conductance background, as 
shown in Figs. |2]c-f. We observe that the ZBCP appears as 
a simple structure in Fig. |2j), while in the remaining spec- 
tra it results to be structured with variations of slope or sec- 
ondary maxima, as in Fig. [2^. The maximum conductance 
ratio G NS {V = 0)/G NN (V » A) is less than 2 for all 
the curves, however G NS {V = 0)/G NN (V » A) ~ 2.2 
for the data in Fig. |2f. In addition to this, the energy width 
of the main zero bias triangular structure is lower than lOmV^ 
in Figs. [2^-d while it results wider, around 40m V, in Figs. 
|2t,f. At a first qualitative analysis, these data appear quite 
puzzling and could be interpreted in term of local, large vari- 
ations of the superconducting energy gap. In the following, 
we will show that the theoretical fittings of all the spectra give 
clear indication of a d-wave symmetry of the superconducting 
order parameter, with consistent values of the inferred ampli- 
tude of the energy gap. 

First of all, let us quantitatively analyze the curves of Figs. 
|2^-d. We were not able to reproduce the conductance spectra 
reported in Fig. |2]a,b by using either the conventional s-wave 
model or the anisotropic one, even by considering small Z val- 
ues, indicative of low barriers. On the other hand, as can be 
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observed in Fig. 1, the s-wave fittings cannot model the struc- 
tured conductances reported in Figs. [3]c,d. The solid lines 
in the figures are the theoretical fittings obtained by consider- 
ing a d-wave symmetry of the order parameter in the modified 
BTK model, Eqs. 12161 A satisfactory agreement is obtained 
by using as fitting parameters the superconducting energy gap 
A, the barrier strength Z, the angle a and a phenomenological 
factor Y Dynes ~ to take into account pair breaking effects and 
finite quasiparticle lifetime 16 . We notice that in the considered 
spectra, both quasiparticle tunneling and Andreev reflection 
processes take place, since intermediate Z values have to be 
used to simulate the barrier strength (0.45 < Z < 0.9). More- 
over, the angle a varies between 0.39 and 0.51, indicating that 
the average transport current mainly flows along an intermedi- 
ate direction between the nodal one (a = 7r/4) and that of the 
maximum amplitude of the energy gap (a = 0). The modified 
d-wave BTK model allows to satisfactorily reproduce the vari- 
ations of slope around ±lmV of the structured ZBCP in Figs. 
|2k,c,d with a light discrepancy in modeling the full height of 
the peak in Fig. [2ji. We show in the next section that a more 
satisfactory fitting for this contact can be obtained by taking 
into account an additional in series intergrain junction. 

We remark that the values of the superconducting energy 
gap, inferred from the theoretical fittings, are all consistent 
and enable us to estimate an average value of the amplitude of 
the order parameter A = (2.8 ± 0.2) meV. This value is sur- 
prisingly low in comparison with the amplitude of the energy 
gap in other cuprate superconductors, however the possibility 
that the presence of the R11O2 magnetic planes can play an 
important role in the complex Ru-1212 system has to be taken 
into account. We notice that the ratio between the smearing 
factor Y Dynes and the superconducting energy gap results al- 
ways less than 20% and it vanishes for the fitting shown in 
Fig. |2p. We consider this fact as an indication of the good 
quality of our point-contact junctions. 



IV. ROLE OF THE INTERGRAIN COUPLING 

To complete our discussion about the spectra measured at 
low temperatures, we now address the analysis of the conduc- 




FIG. 3: (Color online) Intergrain coupling effect in poly crystalline 
samples. The measured voltage Vmeasured is the sum of two terms: 
Vpc, the voltage drops between tip and sample, the N/S PC junc- 
tion, and Vj, the voltage drops between two superconducting grains, 
forming the S/I/S Josephson junction. 



tance curves reported in Figs. |2^,f, with a wider ZBCP. In this 
respect, we observe that, due to the granularity of the com- 
pound, in same cases, an intergrain Josephson junction can be 
formed in series with the Point Contact one, as schematically 
drawn in Fig. [5] This topic has been recently addressed in 
PCAR studies on MgCNi 3 17 and MgB;A 

To provide a quantitative evaluation of the conductance 
spectra, we consider a real configuration in which the Pt-Ir 
tip realizes a PC junction on a single Ru-1212 grain, which, in 
turn, is weakly coupled to another grain, so forming a Joseph- 
son junction. In this case the measured voltage corresponds to 
the sum of two terms: 



V, 



measure' 



d (I) = V PC (I) + Vj(I) 



(8) 



where Vpc and Vj are the voltage drops at the N/S Point Con- 
tact junction and at the S/I/S intergrain Josephson junction, re- 
spectively. This last contribution can be calculated by the Lee 
formula 19 which, in the limit of small capacitance and at low 
temperatures, reduces to the simplified expression 20 : 
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At the same time, for the Point Contact contribution, we use 
again the extended BTK model for a d-wave superconductor. 
The I(V) characteristic is then calculated by inverting Eq. (|8|i 
and the conductance spectrum is given by: 



^ = dV 
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dVj 
~~dT 



(10) 



By applying this simple model we have satisfactory fitted 
the experimental data reported in Figs. |2^,f. Remarkably, 
for both spectra, the best fittings have been obtained by us- 
ing A = 3.0 meV, consistently with the average value ex- 
tracted from the other curves in Figs. [2^-d. We observe that, 
in this model, two more parameters are needed, namely the re- 
sistance Rj and the critical current Ij of the Josephson junc- 
tion. However, the choice of these two parameters is not com- 
pletely arbitrary, since the condition Rj + Rpc — Rnn has 
to be fulfilled, where Rnn is the measured normal resistance 
and the product R jlj necessarily results lower than A 21 . 

In some cases, it has been pointed out that dips in the con- 
ductance spectra can be related to the presence of intergrain 
junctionsiii^, and for sake of completeness, we have applied 
our model also to the spectra of Figs. |2]a-d. We notice that 
for different junctions, the effect of the intergrain coupling 
results more or less evident, depending on ratio Rj/Rpc- 
For the conductances shown in Figs. [2^-c this effect turns 
out to be negligible, however, some improvement of the the- 
oretical fitting is obtained in the case of Fig. |2jl (see dashed 
line). Remarkably, by this last fitting we have found the same 
value of the superconducting energy gap previously inferred, 
A = 2.8meV, with a r D y nes /A ratio less than 3% and 
Rj/Rpc = 0.05. 
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FIG. 4: (Color online) Temperature evolution of the conductance 
spectrum of Fig. |2j: from T = 4.2K to up the critical tempera- 
ture (dots). The solid lines are the theoretical fittings obtained by a 
modified d-wave BTK model with the energy gap as only free pa- 
rameter. 
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FIG. 5: (Color online)Temperature dependence of the superconduct- 
ing energy gap as inferred from the theoretical fittings shown in Fig. 
|4| The solid line is a guide for the eyes. The right hand scale refers 
to the temperature evolution of the measured height of the ZBCP 
normalized to the 4.2K value. 



V. TEMPERATURE DEPENDENCE OF THE 
CONDUCTANCE SPECTRA 

To achieve information on the temperature dependence of 
the superconducting energy gap in the Ru-1212 system, in this 
section we analyze the temperature behavior of the conduc- 
tance spectrum shown in Fig. [2ji. Indeed, this PCAR junction 
resulted to be very stable for temperature variations. 

In Fig. |4]we show the conductance characteristics measured 
in the temperature range A.2K < T < 35K. We firstly no- 
tice that the ZBCP decreases for increasing temperature and 
disappears at about T ~ 30K, that we estimate as the local 
critical temperature T l c of the superconducting Ru-1212 grain 
in contact with the Pt-Ir tip, consistently with the resistivity 
measurements 11 . This fact provides further evidence that the 
ZBCP is a consequence of the superconducting nature of Ru- 
1212 and is not due to spurious effects like inelastic tunneling 
via localized magnetic moments in the barrier region 22 * The 
experimental data for each temperature are then compared to 
the theoretical fittings calculated by using the d-wave modi- 
fied BTK model with a small contribution of Josephson junc- 
tion in series. For all the curves, we fixed the strength of the 
barrier and the angle a to the values obtained at the lowest 
temperature. 

The resulting temperature dependence of the superconduct- 
ing energy gap A(T) is reported in Fig. [5] where vertical 
bars indicate the errors in the gap amplitude evaluation, that 
increase when approaching the critical temperature. Contrar- 
ily to what expected for BCS superconductors, we observe 
that the energy gap, at low temperatures, decreases rapidly for 
increasing temperatures and goes to zero at T l c in a sublin- 
ear way. We notice that the same temperature evolution for 
the superconducting energy gap is found trough the d-wave 
BTK model with or without considering any intergrain junc- 
tion in series; remarkably, in this last case, the superconduct- 
ing energy gap A remains the only varying parameter. A sim- 
ilar temperature dependence has been reported by G. A. Um- 
marino et alQ, however these authors give a larger estimation 
of the maximum gap amplitude. 

From the average value of the superconducting energy gap 
A = 2.8 meV and from the measured local critical tem- 
perature T\ ~ 30K, we obtain a ratio 2&/(k B T l c ) ~ 2 
much smaller than the predicted BCS value and also smaller 
than the values found for high-T c cuprate superconductors 24 . 
Again we speculate that the simultaneous presence of super- 
conducting and magnetic order is an important key for un- 
derstanding the behavior of the Ru-1212 system. Coexis- 
tence of superconductivity and antiferromagnetism is found 
among cuprates, however it is common believe that ferromag- 
netism and superconductivity are mutually excluding orders. 
Recently, it has been found that in conventional Supercon- 
ductor/Ferromagnetic (S/F) structures, proximity effect give 
rise to an oscillatory behavior of the superconducting T c as a 
function of the thickness of the F layej 2 ^ 2 ^. There are con- 
ditions for which a change of sign of the order parameter oc- 
curs, producing the it -junction phenomenon^! In addition to 
this, a dramatic suppression of the amplitude of the order pa- 
rameter is expected for high T c superconductors in close con- 
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tact with a ferromagnetic material and various examples of 
anomalous temperature behavior are found in the literature. 
Gapless superconductivity can be achieved, that can induce a 
sublinear temperature dependence of the superconducting en- 
ergy gap. In the Ru-1212 system, it has been proved that the 
RuC>2 planes are conducting, however these do not develop 
superconductivity at any temperaturaS. By means of differ- 
ent experimental techniques, it has been inferred that a large 
fraction of the charge carriers is not condensed in the super- 
conducting state even at low temperatures 29 . Both findings are 
consistent with a reduced value of the 2A/ (ksT,.) ratio in this 
compound. 

In Fig. |5]we also report (righthand scale) the temperature 
evolution of the height of the ZBCP normalized to its value at 
T = &.2K. It is worth to notice that G NS (V = 0, T), as di- 
rectly measured from the experiments, and A(T), as inferred 
from the theoretical fittings, show a similar scaling with tem- 
perature. This correspondence is easily verified for Z = 
in case of a s-wave superconductor, however it is a quite new 
result since it has been found for intermediate barriers and un- 
conventional symmetry of the superconducting order parame- 
ter. 



VI. MAGNETIC FIELD DEPENDENCE OF THE 
CONDUCTANCE SPECTRA 

As we already observed, one of the most interesting fea- 
tures of the Ru-1212 is the coexistence of the superconduct- 
ing phase and magnetic order. Indeed, from Nuclear Mag- 
netic Resonance (NMR) 30 ' 31 and magnetization 32 measure- 
ments, it has been found that in this compound ruthenium oc- 
curs in a mixed valence state Ru A+ , Ru 5+ with some higher 
Ru 5+ concentration. The RuO-2 planes, from one side, act 
as charge reservoir for the superconducting C11O2 planes, on 
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FIG. 6: (Color online) Magnetic field dependence of the normalized 
dl/dV vs V characteristics at T = 4.2 K from T to 2 T (dots) 
for the spectra of Fig. I2J5. The full lines are the theoretical fittings 
obtained as discuss in the text. In the inset the magnetic field depen- 
dence of the energy gap is reported. 



the other hand, as observed in Muon Spin Rotation (/iSR) 
experiments 8 , they show quite homogeneous ferromagnetic 
order below T c . A weak interaction between the two order pa- 
rameters, ferromagnetism in the R.UO2 planes and supercon- 
ductivity in the C11O2 planes, has been suggested and recently 
several experiments appear to confirm this hypothesis 2 ^. De- 
spite of the huge experimental and theoretical efforts focused 
on the study of the interplay between superconductivity and 
magnetism, to the best of our knowledge no spectroscopic 
studies in magnetic field of the superconducting order param- 
eter in Ru-1212 have been reported in literature so far. 

In Fig. [6] we show the PCAR spectra measured by apply- 
ing an external magnetic field, parallel to the tip, with in- 
tensity H varying from 0T to 2T. The dl/dV vs V curves 
refer to the contact reported in Fig. I2J5. A reduction of the 
ZBCP for increasing magnetic fields is observed, that in first 
approximation can be reproduced by a phenomenological ap- 
proach. Indeed, addressing the problem of the magnetic field 
dependence of the conductance characteristics is non conven- 
tional superconductors, is a quite difficult task and a complete 
treatment of PCAR spectroscopy in magnetic field would re- 
quire the use of an appropriate density of states in calculating 
the BTK expression for the reflection and transmission coef- 
ficients at the N/S interface. Due to the lack of an analytical 
model, Miyoshi et aL 33 presented a two fluid model to repro- 
duce the effect of normal vortex cores in PCAR junctions in 
conventional superconductors, assuming that the contact area 
contains multiple randomly distributed individual junctions 
(non-Sharvin regime 14 ). These authors propose a simplified 
expression for the conductance, written as a sum of normal 
and superconducting channels: 

GNS t jV) = (1 - h)G NN + hG NS {V) 

where h = H/H C2 and H C2 is the critical field. This ap- 
proach, however, cannot be applied to our experiments since 
we deal with policrystalline, unconventional superconductor, 
exhibiting internal magnetic ordering. In this case, the mag- 
netic induction B is not simply proportional to the external 
magnetic field H and as a consequence the density of vortices 
is not linearly related to H. 

An alternative way to perform a theoretical fitting is ob- 
tained by using an additional pair breaking parameter to sim- 
ulate the effect due to the magnetic field^i^. In this case, 
the total broadening effect T is considered as the sum of 
two terms: T = r°y nes + T ext where r^ nes is the intrin- 
sic broadening due to the quasiparticle lifetime, as used in 
the previous fittings, while T ext mimics the pair breaking ef- 
fect due to the external applied magnetic field. The curve at 
H = 0T (see Fig. [2J}) has been fitted by using the d-wave 
modified BTK model with A = 3.0 meV . For increasing 
magnetic fields we keep constant, in the numerical computa- 
tion, the strength of the barrier Z = 0.9, the orientation angle 
a = 0.51 and the intrinsic Y Dynes — 0.7 meV, while varying 
only two parameters: the energy gap A and the magnetic field 
effect Y ext . We observe that the best theoretical fittings (solid 
lines in Fig. [6j satisfactorily reproduce for any field both the 
height and the amplitude of the measured spectra. In the inset, 
we report the magnetic field dependence of the superconduct- 
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FIG. 7: Normalized conductance curves for the contact of Fig. [2^ 
measured at T = 4.2if in magnetic field up to 2.5 T. When the field 
is switched off, the original spectra are recovered. 

ing energy gap (dots) as extracted from the theoretical fittings. 
The amplitude of the energy gap reduces linearly for H up to 
2T and by a linear extrapolation of the data, we find that the 
energy gap disappears at about H ext ~ 30 T, consistently 
with the estimated critical field reported in Refill. 

We have also studied the effect of the magnetic field on the 
conductance characteristics of the junctions showing wider 
ZBCP, that are formed by two junctions in series. In Fig. 
we report the dl/dV vs V curves measured up to 2.5 T for 
the contacts of Fig. [2^. In this case, we observe that the con- 
ductance curves dramatically change with the application of 
the magnetic field. As discussed in the previous section, the 
Josephson current due to the intergrain coupling is immedi- 
ately suppressed by the magnetic field, modifying the spectra 
towards the narrower, non-structured, triangular shape of the 
ZBCP. In addition to this, the oscillatory behavior of the back- 
ground, due to the intergrain coupling disappears in magnetic 
field. Remarkably, for the junctions of both Figs. 16 17 1 the pe- 
culiar features of the spectra together with the normal junction 
resistance, are recovered when the magnetic field is switched 
off, and no hysteresis is found for increasing/decreasing fields. 

VII. CONCLUSIONS 

We have analyzed the PCAR conductance spectra obtained 
in superconducting RuS^GdQ^Os (Ru-1212) policrystalline 



pellets. All the conductance curves at low temperatures show 
a Zero Bias Conductance Peak that decreases for increasing 
temperatures and disappears at the local critical temperature 
T\ ~ 30K of the superconducting grain in contact with the 
Pt-Ir tip. The triangular shape of all the measured spectra has 
been modeled by using a modified BTK model for a d-wave 
symmetry of the superconducting order parameter. This find- 
ing suggests a closer similarity of the Ru-1212 system to the 
high T c cuprate superconductors rather than to the magnetic 
ruthenate Sr2Ru04 compound. However, the remarkably low 
values of the energy gap A = (2.8 ± 0.2)meV and of the ra- 
tio 2 A/fc^Tc ~ 2 indicate major differences between the Ru- 
1212 and the high T c cuprates. We speculate that the presence 
of ferromagnetic order within the superconducting phase re- 
sults in an effective reduction of the energy gap. We have also 
demonstrated that, when dealing with granular samples, inter- 
grain coupling effects can play a predominant role. In some 
cases, an intergrain Josephson junction in series with the point 
contact junction is formed. Taking into account this feature as 
well, all conductance spectra have been properly modeled by 
considering a d-wave symmetry of the order parameter, with 
consistent values of the amplitude of the energy gap. 

By fixing all the fitting parameters to their values at the 
lowest measured temperature, and by varying A, the tempera- 
ture dependence of the energy gap has been extracted from the 
conductance characteristics of a very stable junction. We have 
found that the energy gap exhibits a sub-linear dependence in 
temperature. The magnetic field behavior of the spectra has 
been also studied, showing a linear reduction of the energy 
gap for fields up to 2 T, from which a critical field Hc2 ~ 30 
T is inferred. We have found that both the superconducting 
features and the normal background in the conductance spec- 
tra do not show any hysteresis in magnetic field. These obser- 
vations seem to suggest a weak coupling between the super- 
conducting and magnetic order parameter. 

Our analysis may be helpful for a deeper understanding of 
the mechanisms enabling high temperature superconductivity, 
and its interplay with magnetic order in unconventional super- 
conductors like rutheno-cuprates. 
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